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ABSTRACT 

I describe a method for synthesizing photometric passbands for use with current 
and future X-ray instruments. The method permits the standardisation of X-ray 
passbands and thus X-ray photometry between different instruments and missions. 
The method is illustrated by synthesizing a passband in the XMM-Newton EPIC 
pn which is similar to the ROSAT PSPC 0.5-2 keV band. 

Key words: methods: data analysis - methods: observational - techniques: pho¬ 
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1 INTRODUCTION 

Photometry, the measurement of source brightness in 
a given wavelength range, is one of the most funda¬ 
mental techniques in astronomy. It is practised through¬ 
out the electromagnetic spectrum, and X-ray astron¬ 
omy is no exception. In optical astronomy, photomet¬ 
ric passbands have been standardised between observa¬ 
tories since the introduction of the John son and Mor - 
gan system (IJohnson fc Morgan. 195~3l ); see iBessel (20051) 
for a review. The benefits of a well-defined set of pho¬ 
tometric X-ray bands have been recognised before, (e.g. 
IGrimm et al., 20091) but despite more than 50 years of X- 
ray astronomy, there is as yet no system of photometric 
passbands that applies to more than a single X-ray obser¬ 
vatory. In this letter I will describe the concept of synthe¬ 
sized X-ray passbands that can be applied to existing and 
future X-ray astronomy instruments, and so permit the es¬ 
tablishment of standard, multi-mission X-ray photometric 
passbands. As an example, I will show that the ROSAT 
PSPC 0.5-2 keV passband can be synthesized with the 
XMM-Newton EPIC pn camera, so that X-ray photome¬ 
try gathered with ROSAT and XMM-Newton can be com¬ 
pared directly. 


2 DEFINITIONS 

I will adopt a definition of photometry as the measurement 
of source brightness through a well defined passband rel¬ 
ative to some standard reference. The standard reference 
could be an astronomical object such as the star Vega. Al¬ 
ternatively it could be a theoretical reference, such as the 
constant flu x density spectru m used as a reference in the 
AB system dOke. 1965lf 19741 ). provided that the through¬ 
put of the instrument performing the photometric mea¬ 
surements is sufficiently well calibrated. Almost all mea¬ 
surements of source brightness in X-ray astronomy may be 


considered to conform to such a definition: the curves of ef¬ 
fective area against energy commonly employed in X-ray 
astronomy represent well defined passbands which have 
already been calibrated in physical units. However, where 
X-ray photometry is disadvantaged with respect to e.g. 
optical photometry, is that each instrument on each X- 
ray observatory has a different effective area curve. To 
mitigate this problem, X-ray astron omers have developed 
tools such as PlMMif] dMukai. 19931 ) which allow a mea¬ 
surement with one instrument to be translated into a pre¬ 
dicted count rate for a different instrument, given a model 
for the X-ray spectrum of the source. However, if the spec¬ 
trum of the source is not known or is only poorly deter¬ 
mined, then the predicted count rate will not be accurate. 

In a photon counting detector, as usually employed in 
X-ray astronomy instrumentation, the count rate C can be 
related to the spectrum of a source and the effective area 
of the instrument as follows: 

o-r^ m 

" ^min 

where v is frequency, h is Planck’s constant, fv(v) is the 
flux per unit frequency interval of the source and A(u) 
is the effective area of the instrument. For an imaging 
system, it is assumed that A(y) is the overall effective area 
of the instrument and telescope combined. The limits Vmin 
and Vmax are the minimum and maximum frequencies at 
which A(v) is significant. 

If we define 


K = max 


My) 

hv 


and 

R{*) 


A{u) 

Khv 


1 https://heasarc.gsfc.nasa.gov/docs/software/tools/pimms.html 
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so that R{u) is the normalised respose curve of the system, 
then 

^ = f f„{v)R{v)dv (2) 

J Vmin 

The fundamental quantity measured in astronomical 
photometry is the quantity on the right hand side of Equa¬ 
tion [2] the integral over the passband of the flux of the 
source multiplied by the response of the system. The quan¬ 
tity on the left hand side of Equation [2] is the measurement 
itself, the count rate divided by the peak effective area of 
the system. In optical astronomy the measurement is usu¬ 
ally expressed as a magnitude m, 

r^max 

m = — 2.51og 10 / f v (v)R(v)dv + z 

J Vmin 

where z is the zeropoint, defined such that 
a reference standar d has zero magnitude (e.g. 
ICousins fc Jones. 19761 1. A Vega-based magnitude system 
would not be appropriate in the X-ray regime, because 
Vega has yet to be convincingly detected as an X-ray 
source ((Ayres, 20Qg| ; [Pease. Drake fe Kashvap, 2006fo , but 
the AB magnitude system would be well suited to X-ray 
astronomy. Nonetheless, the adoption of magnitudes as a 
convenient unit is not necessary for the establishment of 
a system of photometric passbands in X-ray astronomy. 
What is essential is the ability to perform the measure¬ 
ment described by Equation [2] with different instruments 
on different observatories with similar response functions, 
or passbands, R(v). 


3 SYNTHESIZING PHOTOMETRIC 

PASSBANDS 

The synthesis of a passband is only possible for instru¬ 
ments which discriminate the energies of the incoming 
photons. In X-ray astronomy, data are usually telemetered 
to Earth as ‘event lists’, in which the individual photons 
are listed with their properties such as arrival time, posi¬ 
tion on the detector and energy channel. The energy res¬ 
olution of non-dispersive Charge-Coupled Device (CCD) 
based X-ray instruments is limited, usually 1 0 < R < 50 
(e.g. iTurner et al., 200ll ; IStriider et al., 200 jl . so the re¬ 
lationship between energy channel and (real) photon en¬ 
ergy is described by the redistribution matrix. The redis¬ 
tribution matrix contains, for a sequence of discrete en¬ 
ergy ranges, the probability of a photon within that energy 
range being recorded in each of the energy channels. I will 
refer to the product of the redistribution matrix and the 
effective area as a function of energy (i.e. for each discrete 
energy range, the channel probability distribution multi¬ 
plied by the effective area in that energy range) as the 
response matri>0. For each energy channel i, the response 
matrix contains the effective area for that channel as a 
function of energy, and therefore frequency Ai(v). 

X-ray astronomers choose passbands for their instru¬ 
ments by selecting ranges of energy channels which corre¬ 
spond to their desired energy ranges; the Chandra ACIS 


2 It should however be noted that, somewhat confusingly, it is 
quite common in X-ray astronomy to use the term ‘response 
matrix’ interchangably with ‘redistribution matrix’. 


bands advocated by [Grimm et al. ('20091'). and the bands 
chosen for the 2XMM catalogues ( Watson et al., 20091 1 are 
examples of this approach, and I will refer to these bands 
as natural instrumental passbands. For natural passbands, 
the count rate in the band C is related to the count rates 
in the individual energy channels Ci by 


where i m i„ and imax are the minumum and maximum 
channels of the passband. The effective area of the band 
A(u) is related to the effective areas A;( v) of the individual 
channels by 

A W = M v ) 


While natural passbands can be chosen to cover similar 
energy ranges with different instruments on different mis¬ 
sions, there is little scope in this approach to control (and 
thus standardize) the shapes R{v) of the passbands. 

To synthesize a passband, it must be possible to con¬ 
trol the shape as well as the energy range. This can be 
accomplished by assigning weights wi to the channels that 
form the band when calculating the count rate within the 
band. Thus for a synthesized band: 

C = V w-iCi 


and 

A(u) = ^ w-iAi{v ) 


The weights can be chosen so that the synthesized band 
has a response R{v) as close as possible to the response 
desired. The weights Wi therefore act as a synthetic filter, 
controlling the shape of the passband. The capability to 
compose and apply the synthetic filter after the data have 
been taken is a distinct advantage of X-ray instrumenta¬ 
tion compared to that employed for photometry at longer 
wavelengths. 


4 PRACTICAL CONSIDERATIONS 

In order to synthesize a passband with an instrument, the 
instrument must be sensitive over a similar or larger fre¬ 
quency range than that covered by the passband. The de¬ 
gree to which the shape of a passband can be controlled by 
the choice of the weights Wi depends on the width of the 
channel distribution at a given energy in the response ma¬ 
trix (i.e. the energy resolution of the instrument) and to a 
lesser degree on the shape of this distribution. The intro¬ 
duction of the weights Wi implies a reduction in the signal 
to noise of the count rate with respect to that in a natural 
passband. In general, X-ray instruments have improved in 
both spectral resolution and collecting area with time, and 
this trend is expected to continue with the future Euro¬ 
pean Space Agency mission Athena (iNandra et al., 20131 '). 
Thus it is more practical to use present-day or future in¬ 
struments to synthesize the natural bands of earlier X-ray 
astronomy instruments than the other way round. 

In general, the more closely the response curves R(y) 
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of two passbands match, the more similar the photome¬ 
try will be between them. Therefore the best photometric 
performance from a synthesized band will be obtained by 
choosing the weights Wi to match as closely as possible 
the desired R(u). On the other hand, the smoother, and 
more uniform the Wi, the higher the signal to noise will 
be. Therefore, in choosing the weights w,, there might be 
some trade-off between obtaining the desired shape R(v) 
and maximising the signal to noise ratio, particularly if 
the best match to R(v) corresponds to large fluctuations 
in Wi between adjacent channels. 

Another consideration is that the count rate in a syn¬ 
thesized band will no longer be Poisson distributed as it 
is in natural bands; instead it will be distributed as the 
sum of the weighted Poisson distributions, and error anal¬ 
ysis is more involved. However, the subtraction of back¬ 
ground also causes the count rates of sources to deviate 
from Poisson statistics, and this is routinely performed in 
X-ray astronomy already, so the importance of this con¬ 
sideration should not be exaggerated. Nonetheless, the re¬ 
duction of signal to noise and the increased complexity of 
error analysis are disadvantages associated with synthe¬ 
sized passbands, to be traded off against the advantages 
of standardising the passbands between missions. 


5 AN EXAMPLE SYNTHESIZED PASSBAND 

As an example, I will now show that the 0.5-2 keV 
band of the ROSAT Position Sensitive Proportional 
Counter can be synthesized with the XMM-Newton Eu- 
ropean Ph oton Imag ing Camera (EPIC) pn instrument 
dStriider et al., 20011 ). To date the only imaging all-sky 
survey carried out in X-rays used the ROSAT PSPC, so 
this instrument has provided p hotometric measur ements 
of a large body of X-ray sources dVoges et ah, 19991 '). Many 
of these sources have been observed subsequently with 
XMM-Newton. It would be useful to be able to compare 
directly the photometric measurements from ROSAT and 
XMM-Newton, for example to study the long term vari¬ 
ability of X-ray sources. Following the logic of Section |4j 
it is more practical to synthesize the ROSAT PSPC pass- 
band with the XMM-Newton EPIC pn than the other way 
around, because EPIC pn has better energy resolution, 
and a larger effective area than the PSPC. 

The ROSAT PSPC 0.5-2 keV band is formed from 
PSPC energy channel^] 52 - 201. The normalised response 
R(u) of this band is shown as the grey curve in Fig. [T] The 
distribution of weights Wi to be applied to the EPIC-pn 
channels was obtained by iteratively adjusting the weights 
so as to minimize the sum of the squares of the difference 
between the PSPC and synthesized response R(v) at each 
point of the R(v) curve over the frequency range of the 
synthetic band, i.e. a least-squared fit. No smoothing crit- 
era were applied during the fitting of the weights. Fig. [2] 
shows the channel weights Wi corresponding to the syn¬ 
thesized passband. The distribution of Wi is quite smooth 

3 These are usually called pulse-height invariant or PI channels, 
because they correspond to the pulse height of an event after 
it has been calibrated for temporal and spatial gain variations 
onto a standard energy scale 



Figure 1. Passband of the ROSAT PSPC 0.5 - 2 keV band 
(solid curve) and an XMM-Newton EPIC-pn passband (dashed 
curve) which has been synthesized to have approximately the 
same shape. 


with fluctuations between adjacent channels rarely exceed¬ 
ing 2 per cent, so smoothing the Wi would have little im¬ 
pact on the signal to noise ratio of photometry through 
the synthesized band. The synthesized R(v) is shown as 
the black, dashed curve in Fig. |T] The match between the 
ROSAT PSPC and EPIC pn synthesized bands is good, 
with R(v) differing by less than 0.01 over most of the band, 
and by no more than 0.05 at any point. However, the weak 
tail in the synthesized passband which stretches to higher 
frequency is irreducible, because it corresponds to the low 
energy tai0 of the CCD response. 

To investigate the performance of the synthesized 
band in reproducing PSPC photometry, I have similated 
a number of spectral models and folded them through the 
responses of the PSPC, EPIC-pn_and the sy nthesized 0.5- 
2 keV band using XSPEC11 (jArnaud. 19961 ). For the nor¬ 
mal EPIC-pn 0.5-2 keV response, I simply selected the pn 
channels that have nominal energies within the range 0.5- 
2 keV. The spectral models were chosen to span the range 
of spectral shapes typically exhibited by sources in the 0.5- 
2 keV band, and comprise power laws, with photon index 
F =1, 2 and 3, a power law with F = 2, absorbed by cold 
gas with a column density of Nh = 10 22 cm~ 2 , optically 

4 i.e. redistribution of photons of a given energy to channels of 
lower energy 
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Table 1. Simulated 0.5-2 keV photometry in the PSPC, EPIC-pn, and the synthesized PSPC band, for a number of spectral models. 
The second, third and fourth columns give the 0.5—2 keV flux that would be observed if the given spectrum is observed using the 
PSPC, EPIC-pn and Synthesized PSPC bandpass respectively, and the count rate is translated to a flux using a standard energy 
conversion factor. In all cases the spectral model was normalised such that the true 0.5-2 keV flux is l.OOxlO -14 erg cm~ 2 s' 1 . 
The last two columns give the ratios of the flux that would be observed with EPIC-pn and the synthesized band respectively, to the 
flux that would be observed with the PSPC. 

Model — Simulated Flux (10 -14 erg cm -2 s _1 ) — EPIC-pn / Synthesized / 

PSPC EPIC-pn Synthesized PSPC PSPC 


power law T = 1 

0.96 

power law F = 2 

1.02 

power law F = 3 

1.04 

power law F = 2, Nh = 10 22 cm -2 

0.68 

mekal kT = 0.3 keV 

1.11 

mekal kT = 0.7 keV 

1.32 

mekal kT = 1.5 keV 

1.16 

black body kT = 0.3 keV 

1.09 



Figure 2. Weighting function applied to the EPIC pn channels 
to produce the synthesized band shown in Fig. |T] 

thin thermal plasmas (mekal in xspec) with kT =0.3, 0.7 
and 1.5 keV, and a black body with kT =0.3 keV. Each 
of the spectral models was normalised so as to have the 
same true 0.5-2 keV flux of 10 -14 erg cm -2 s _1 , and the 
count-rates were simulated with zero statistical error, so 
that any differences in the photometry from the 3 pass- 
bands are systematic differences related to the shapes of 
of the passbands. 

In X-ray astronomy, photometry is normally derived 
from imaging surveys by measuring the count rate in a 
given band, and dividing this by an “energy conversion 
factor” (ECF) which is specific to the instrument and band 


0.96 

0.96 

1.00 

1.00 

1.02 

1.02 

1.00 

1.00 

1.04 

1.03 

1.00 

0.99 

0.84 

0.71 

1.25 

1.04 

1.11 

1.11 

1.00 

0.99 

1.12 

1.31 

0.85 

0.99 

1.03 

1.15 

0.89 

0.99 

1.01 

1.08 

0.93 

0.99 


to obtain a flux in physical units. We have followed this 
process in simulating the photometry that would be ob¬ 
tained through the different passbands. The ECF itself is 
typically derived by folding a simulated source of known 
flux and spectral shape through the response of the instru¬ 
ment following the same procedure as we are employing to 
obtain the count rates for our simulated sources. To gen¬ 
erate ECFs for the three passbands, we adopt the spectral 
model used for the 3XMM catalogue, a power law of pho¬ 
ton index T = 1.7 absorbed by cold ga s with a c olumn 
density Nh = 3 x 10 20 cm -2 llRosen et al.. 20~i~5h . The 
ECFs for the PSPC, EPIC-pn and the synthesized band 
are 7.65 x 10 10 , 6.67 x 10 11 and 4.75 x 10 11 ergs -1 cm 2 
respectively. 

The results of the simulated photometry are given 
in Table [l] It can be seen that the flux derived from the 
PSPC differs significantly for several spectral models from 
the flux derived from photometry in the natural EPIC-pn 
0.5-2 keV band. In the worst case of the heavily absorbed 
spectrum, the flux differs by 25 per cent. In contrast, the 
flux derived from the synthesized band differs from the 
PSPC flux by less than 1 per cent for all spectral models 
except for the strongly absorbed spectrum. In that case 
the difference is only 4 per cent, and is due to the weak 
tail towards higher energy in the synthesized response, 
where the absorbed spectrum is rising rapidly. 


6 CONCLUSIONS 

I have outlined a method for synthesizing photometric 
passbands in X-ray astronomy. It offers the potential to 
standardize photometry from different X-ray astronomy 
instruments and missions. I have limited this letter to a 
conceptual description of the approach, and provided an 
example to illustrate how it might be used. Widespread, 
practical application of the method will require consider¬ 
ably more work. Software tools will be required to opti¬ 
mize the shapes of the synthetic filters w ; used to shape 
the passbands, and software tools will be needed to com¬ 
pute photometry for real astronomical sources in the syn¬ 
thetic passbands. Not least, for this method to facilitate 
the large-scale provision of photometry in standardized X- 
ray passbands, there will need to be concordance within 
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the X-ray astronomy community, at least amongst the cu¬ 
rators of archives and the constructors of astronomical 
catalogues, as to what to choose as the minimum set of 
standard passbands. 
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